On Brownian flights 



^ ; Athanasios BATAKIS*, Pierre LEVITZtand Michel ZINSMEISTER* 

O 

(N 

5h , 

^ Abstract : Let K be a compact subset of PJ\ We choose at random with uniform law 

a point at distances of K and start a Brownian motion (BM) from this point. We study 
the probability that this BM hits K for the first time at a distance > r from the starting 
point. 
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1 Introduction and motivation. 

> 

Porous materials, concentrated colloidal suspensions or physiological organs such that 
£S) ■ lung or kidney are systems developping large specific surfaces with a rich variety of 

■^j- ■ shapes that influence the diffusive dynamics of Brownian particles. A typical example 

is the diffusion of water molecules in a colloidal suspension. NMR relaxation allows to 
measure the statistics of the flights of these molecules over long colloidal shapes such 
as proteins or DNA chains. It is thus tempting to rely these statistics to the geometry 
of the molecules, the goal being to probe shapes using this method. An ideal (but far- 
reaching) objective would be to make up a DNA-test for example using NMR relaxation. 
This program has been developped in |GKL + 06j where various kind of simulations or 
experiments have shown remarkable commun properties. 

All the simulations measure the statistics of the same random phenomenon: an irregular 
curve or surface is implemented, consisting of a union of a large (but finite) number of 
equal affine pieces. Such a piece is chosen at random with uniform distribution and a 
random walker is started at some small but fixed distance from this piece inside the 
complement of the surface. One is interested in the law of the variable X = the length 
of the flight, i.e. the distance between the starting point and the first hitting point on 
the surface of the random walker. Whatever shape the surface shows, the experiment 
shows the same behaviour 

P (X > r) r d *- d - 2 (1) 
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where d e is the dimension of the ambiant space and d is the Minkowski dimension of 
the surface. 

A first mathematical explanation of this behavior was outlined in |GKL+06| . The mam 
purpose of the present paper is to give first a rigorous statement of this result and to 
prove it with minimal assumptions so that all the cases of the simulations are covered. 
This will be the content of the second paragraph. 

The third section concerns an alternative approach of the result in a special but impor- 
tant case: a quasiconformal perturbation of the line in 2D. The paper being dedicated 
to Fred Gehring is one of the reasons of this section, but not the only one: indeed this 
alternative proof is neat and instructive. 

There is a case of particular importance for physics and particularly polymer physics. 
It is the case of a curve being a self-avoiding walk, since this has been shown to be a 
good model for polymers. In 2D there are two evidences that the result should remain 
in this case : first intensive simulations performed by P.Levitz and secondly some com- 
putations by Duplantier using conformal field theory. We present in the last section 
these simulations and sho how, despite the fact that this case is not covered by results 
of section 3, one can use results about SLE to prove the result in this situation. 
We would like to conclude this introduction with two remarks: 

1) The result may be surprising since it involves only the Minkowski dimension while 
the experiment suggests that harmonic measure is involved (and it is!) and consequently 
that the result should depend on some multifractal property of the harmonic measure. 
This is not the case because of the law of the choice of the starting point. A completely 
different behavior would occur if instead of uniform law we had chosen harmonic mea- 
sure and this case is extremely interesting since essentially it models the second flight: 
this case will be considered in a forecoming paper. 

2) The paper |GKL+06| has been written while the second author was hosted by the 
laboratory of physics of condensed matter at Ecole Polytechnique whose members he 
thanks for their warm hospitality. This paper is an example of a successfull pluridis- 
ciplinar research and we would like to emphasize the fact that firstly problems coming 
from physics (especially polymer physics) are extremely rich mathematically speaking 
and that, secondly, modern Function Theory as we herited from great mathematicians 
as Fred Gehring is an extremely efficient tool to attack these challenging problems. 
Our goal will be achieved if the reader gets convinced of this last statement after having 
read this paper. 

2 Geometric considerations 

For convenience all domains in the next two sections of this paper will be assumed to 
have compact boundary. An estimate like cannot be true for every domain. Some 
geometric conditions are needed : one of them is that Minkowski dimension of the 
boundary exists, i.e. that d = lim^o lo ^^ t e exists, where N £ is the minimal number of 
cubes of size e needed to cover the boundary. 



As we shall see later Whitney decomposition is a central tool in our proof; a second 
condition we have to impose is that the number of cubes of the Whitney decomposition 
at distance r must also be comparable to r~ d . A sufficient condition for this is given by 
the first property of NTA domains ( |JK82j ) called the "corkscrew" condition: 

Definition 2.1 We say that a domain Q satisfies the "corkscrew" condition if there 
exists ro > and a constant c > such that for all r < ro and any x G dQ there exists 
y G Q such that cr < dist(x, y) < r and dist(y, dQ) > cr. 

Proposition 2.1 Under the corkscrew condition the number of Whitney cubes with size 
that intersect the level surface T r = {x G Q ; dist(a;, dQ) = r} is comparable to the 
minimal of cubes of size r needed to cover the boundary. 

The straightforward proof is left to the reader. 

3 The case of open sets in R n , n > 3 

Although our approach can be easily adapted to open sets in the plane, we present it for 
M. n , n > 3 for two reasons: First, Green function and related formulas being different 
we try to avoid writing everything twice. Secondly, a different approach is proposed in 
section H] for open sets in the complex plane, using the quasi- conformal theory. 

Notation 3.1 Given an open set Q, a point 16O and a set A G Cl we denote P(x 
A) the probability that Brownian motion started at x touches A before leaving Q. The 
Green function of Q will be denoted Gq and let G n be the Green function of M n . The 
ball of center x and radius r is denoted by M(x, r) and the distance from x to dQ is 
denoted by d x . We say that two quantities A, B are "comparable" (we denote A ~ B) 
if the exists a constant c such that ^A < B < cA. The harmonic measure of a set F 
(usually a subset of dQ) at a point x G Q is denoted by loq(x, F) or lo(x, F, Q). 

3.1 Green estimates 

Proposition 3.1 Let Q be a domain in W 1 , n > 2 and x, y G Q. There exists a universal 
constant C depending only on n such that for all I < | 




Proof First, notice that by the maximum principle, Gw x ^ x ) ( x > z ) — Gn( x , z ) — G n (x, z) 
for all z G B(x, d x /2). It is then easy to check that there exists a constant c = c(n) such 
that c~ 1 G n (x, z) < Gn(x,z) < G n (x,z) for all z G B(x, £d x ). 



Therefore the function that assigns s i— > (£d y ) n 2 G n (s, y) is harmonic in fi\B(?/, £d y ), 
tends to at dQ and takes values between c _1 and 1 on dM(y,£d y ). The proba- 
bility P^x c — >q M(y,£d y )^ is then comparable to (£d y ) n ~ 2 G n (x,y). In a similar way 

^(y M(x, £d x )^j is equivalent to (£d x ) n ~ 2 G n (y, x) and the proof is complete since G n 
is symmetric. • 

We suppose from now on that the domain Q is uniformly "fat" , i.e. that there exists 
a constant c > such that for any i£ll and r < 1, we have 




In particular, this implies that there exists a uniform lower bound L > of the proba- 
bility that Brownian motion started at x hits the boundary of Q before leaving M(x, 2d x ) 
(see |Anc86] . lemma 5). This notion (also called "uniform capacity density condition") 
has previously been introduced in various contexts, cf. [JW88] . [Anc86] . [HKM93] . 



3.2 Main Results 



We consider a Whitney decomposition of Q in dyadic cubes Q satisfying ci|Q| < 
d(Q,dfl) < C2IQI, where c\ < 1 < C2 are positive constants (powers of 2, for conve- 
nience) depending on n. For t > we note Q t the subcollection of cubes of the Whitney 
decomposition that intersect the level surface r t = {x G Q ; d x = t}. 

Theorem 3.2 Take e < r, fix a cube Q r G Q r and consider, for every Q G Q e , its 
center xq G Q. Then p{%q Qr) is equivalent to (^) n 2 . 

QeQ £ 

Proof According to proposition 13.11 we have 

i E P { x Qr 0) < {^Y~ 2 £ P(x Q «- n Qr) <CJ2 P ( x Qr q), 

where xq t is the center of the cube Q r . We now show that p(xq t Qj is 

QeQe 

equivalent to the harmonic measure un(xQ r , dQ) of dQ at xq t (in Q) which equals 1. 
For this purpose we will use hypothesis (TS]) together with an easy control of multiple 
coverings. 

Take any Q G Q e and consider the cube 3c2<5 of same center but of 3c2 times the 
sidelength of Q (02 being the constant of the Whitney decomposition). By condition (T5]) 
the probability for Brownian motion started anywhere in Q to exit Q before exiting 3c2<5 
is bounded below by a positive constant c. Hence, the harmonic measure ujq(xq t , 3c 2 Qn 



d£l) of 3c 2 Q fl <9f2 at xq t (in Vt) is greater than cP \xq t Qj ■ Summing over all cubes 
Q E Q e we get 

C P ( X Qr ^S] Q)< Yl ^i^Qr^^Q) 

On the other hand, every x G <9f2 can only belong to a finite number of cubes 
3c2<3; this proves that ^Xg r , 3c2<3^ is comparable to the harmonic measure of 

Q&Qe 

the boundary dVt of Vt and therefore P ( x Qr <5 1 has an upper bound depending 

only on n and on condition's (T5]) constant. The lower bound is trivial by the "fatness" 
condition. • 



Theorem 3.3 Choose Q at random with uniform law in Q £ . The probability for a 
Brownian motion started at any point x of Q to hit T r before exiting Q is comparable to 
#a /r\»-2 

#Q £ \e) ' 

Proof First observe that by the Harnack principle we can choose x = xq the center of 
the cube Q. Since the cube Q is arbitrarily chosen there are #<2 e possible choices. 

We consider now the cubes of Q r and we define S r as the part of the boundary of 
^QeQrQ separating the set T r = {x <E Q such that d x = r} and dQ. We say that a cube 
in Q r has a seashore if part of its boundary is also part of S r , and , in this case, this 
access to the "sea" is at least a square whose diameter is > to a constant depending 
only on n times the size of the cube. We then consider the open set U consisting in the 

union of the components of Vt \ [^J Q. We denote by S r the boundary of this open 

QeQr 

set: we are interested in the probability that Brownian motion started at any point x 
of Q hits S r before exiting Q. Denote by V the component of U containing Q. Let O 
be a cube in Q r having a seashore. Each one of its sides touching S r contains a dyadic 
square R of ci/8 times the size of O, c\ being the constant in Whitney decomposition, 
such that R C S r . Let Ol be a cube of the same center as O but 1 + Ci/2 its size. We 
consider the dyadique cube R' contained in Ol of size ci/8 times the size of O vertically 
above R and at distance ci/8 from O . 

We use the Boundary Harnack Principle to prove that the probability that Brownian 
motion started at Xq leave V through R is comparable to p{xq R'\ 

By "adapted cylinder" to a graph of a Lipschitz function we understand a vertical 
revolution cylinder of finite height centered on the graph. Let us remind the Boundary 
Harnack Principle : Let u and v be positive harmonic functions on a Lipschitz domain 
vanishing on the graph between the adapted cylinder (to a graph-component of the 
boundary) C and the "sub" -adapted cylinder C of the same center and revolution axis 



but of i times the size, I < 1. If C is the "middle" cylinder of the same center and 
revolution axis but of times the size of C. Then for all x G dC fl V 



v(x) v(P) 



where P is the intersection point of the revolution axis of the cylinder C and of its 
boundary. The multiplicative constants in the equivalence relation depend on the ratio 
(revolution radius): (height) of C , on £, on the local Lipschitz norm of the boundary and 
on the dimension of the space n, see |Anc78] . 

Remark that Ol only touches the neighbouring cubes of R. Clearly, Ol H V is a 
Lipschitz domain (its boundary is composed of a finite union of squares). We can find a 
finite number of adapted cylinders Ci that do not touch the cubes R' such that the "sub"- 
adapted cylinders Ci of half their size cover the boundary of OR fl dV . Furthermore, 
we can ask these cylinders to touch the boundary of the cube 0/ of the same center as 
O but of 1 + ci/100 the size (see figured]). Note that, since there is a finite number of 
configurations of the neighborhood of O in the Whitney decomposition, the number of 
adapted cylinders needed to this covering is bounded by a uniform constant k depending 
only on n. We consider the harmonic functions u(x) = u(x, R', V \{J " we t" sides of o 
and Vi(x) = lu(x, Q fl 80l, V) in the domain V = V \ [Ju we ^> g^es Q f R' ■ Clearly, 



We can apply the Boundary Harnack Principle to every one of these cylinders Ci for 
the functions u and Vi. We get that for all i and all x G dCi fl V 



where Pi is the intersection point of the revolution axis of the cylinder and of it's bound- 
ary. The boundary of our domain is the inner boundary of the cubes of the Whitney de- 
composition; hence we may restrain ourselves to a finite number of configurations of the 
"adapted cylinders" (up to a contractions-dilatations that do not affect constants) and 
therefore the multiplicative constants in the equivalence relation will be finite in number 
and so uniformly bounded away from and infinity. Since the distance of dist(Pj, R') is 
equivalent to the distance dist(Pj, OR) we can prove by standard arguments of harmonic 
analysis that u(P.j) is bounded below by a constant depending only on dimension n. 
Using the Harnack principle and the fact that all adapted cylinders intersect dRi we 
get the existence of a constant c depending only on n such that Vi(Pk) < cVi(Pj) for all 
i,j,k. 

After summing over i, taking in account equations and (jlj) and using standard 
Harnack inequalities we get 



(x, 80 n dV, V) < v i( x ) < KUJ (x, dO n dV, V). 



(3) 



Vi(x) Vi(Pi) 

u(x) u(Pi) ' 



(4) 



to(x, dR n dV', V) < u(x, R', V) < k'w(x, dR n dV\ V) 



o 
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V 



Figure 1: A configuration in IR 2 . 

for all x G dRi U UidCi R V. By the maximum principle we thus obtain 
u){x, OR n ST/', V) < uj(x, R', V) < k'u(x, or n <9V, V") 
for all x e Q e Q e - 

Theorem 13.21 applies to the domain V where Q r is replaced by the collection R' 
coming from all "wet" sides of all cubes with a seashore. We obtain that 

£ p (x Q ^v> R') ~ Q ^ (x w , on n V) 
Qe Q £ , QcV 

a quantity bounded from below by some a > by the fatness condition. The proof 
is completed by suming over all R' and V's and by noticing that up to a multiplica- 
tive factor depending only on the dimension, the number of cubes with a seashore is 
comparable to the number of cubes in Q r . • 



Theorem 3.4 Choose Q at random with uniform law in Q £ . The probability for a 
Brownian motion started at any point x of Q to exit Q at distance greater than r from 

the starting point is comparable to — ^ ( —\ 



Proof It depends on a comparison between the probability of "cruising along the coast" 
dVt and the probability to move at distance r before coming back to the coast. The sec- 



ond is comparable to (^ _ J according theorem 13 . 31 while the first is exponentially 

small. 

Take s > and xq G Q G Q £ . Consider the annuli centered at x of inner radii is and 
outer radii (£ + l)s where £ = 0, [£] . Brownian motion started at x and moving at 
distance r from x before exiting Q must go through all these annuli. The probability of 
going through such an annulus while staying at distance at most | from the boundary 
is bounded by a po G (0, 1) by the "fatness" hypothesis. To see this take any point 
y in the middle of the annulus (i.e. at distance ^-s from x) and consider the ball of 
center y and radius |. If d y < |, the probability to exit the ball without touching dQ 
is uniformly bounded away from 1 by the "fatness" hypothesis. This probability being 
greater than the probability of going through the annulus we have the statement. By the 
independence of the "crossing annulli" events we get that the probability that Brownian 

motion goes through all the annulli is smaller that Pq . 

Let us now prove the following statement: "there exist < Ci,C2 < 1 positive 
constants depending only on dimension and on the constant L of the fatness condition 
such that for any x G Q there exist disjoint sets K\,Ki C dQ H M(x, 2d x ) verifying 
dist(_Ki, K-i) > c\d x and u>q(x, Kj) > C2 for i — 1,2. Once more this is a consequence 
of the "fatness" property. Cut the sphere dM(x, 2d x ) in small equal normal polygons 
(spherical triangles in R 3 ) and consider the intersections Lj of the cones of summit 
x and basis these polygons with the set ~B(x,2d x ) \ M(x,d x ) (see figure [2]). Clearly, 
we can choose the polygons small enough (independently of x and d x ) to have that 
u(x, Li, M(x, 2d x ) \ < L/2 n , where L is the constant in the fatness condition. It 
is now clear that u(x,Li R dQ, M(x, 2d x ) fl Q) < L/100 and therefore there are two 
non-neighboring Li's having harmonic measure greater than L/2 n #Lj, which proves the 
statement. 

The previous statement implies that once we reached distance r from the boundary 
of the domain the Brownian motion will revisit the boundary at distance comparable 
to r from the starting point with probability greater than c<i- Putting together all the 
above we get the theorem. • 

As a corollary we get the following. 

Theorem 3.5 If Q satisfies 

1. The corkscrew condition 

2. The fatness condition 

3. and if dQ has a Minkowski content, 

then for every rj > there exists a constant c v ^ n > such that for all r > e > if we 
choose Q at random with uniform law in Q e , the probability that a Brownian motion 



Figure 2: Note that d(x) is the distance of X from the boundary K = dVt. 

started at any point x of Q hits for the first time dfl at distance greater than r from the 
starting point P (X > r) verifies 



4 An alternative approach in the quasicircle pertur- 
bative case. 

We present here a simple 2D case for which the proof of the main result is particularly 
simple using conformal mapping which preserves Brownian trajectories. The curve we 
consider will be a quasiconformal perturbation of the line. More precisely we consider 
a domain Q = where ip : — > C is holomorphic and such that 



It is known that under this hypothesis ip is injective and has a quasiconformal extention 
to the whole plane. In particular T = (p(M) is a quasicircle close to a line. We will 
moreover assume that T has a Minkowski dimension which we denote by d. By Koebe 
theorem, the quantity y\(p'(x+iy) \ is uniformly comparable to the distance from ip(x+iy) 
to T. For a > we then define L a = {x + iy G IR+ ; y\(p'(x + iy)\ = a} and (p(L a ) will 
serve as a substitute for the level set {( G fl ; d((, T) = a}. 




X 



cp"(x + iy) 
sup y\ — 



< 1/2. 



(5) 



Lemma 4.1 If ' tp satisfies then L a is the graph of a Lipschitz function f a : R — > R. 

Proof We apply the implicit function theorem to the function 

F(x, y) = y<p'(x + iy)<p'(x + iy) — a 2 . 

The computation gives 

OF dF 

— = 2y\tp'{x + iy)\ 2 U^{x + iy)), — = \<p'(x + u/)| 2 (l + 2yT{^{x + iy)) 

where 

ip'(z) 

The result follows because ^ > and by (jSJ). 

We now consider a portion of the curve with diameter 1 and < e < r. We divide 
the portion of Q between T and L £ into pieces of diameter ~ e. The preimages of 
these pieces are rough squares of sidelength f E (xj) where Xj is the left-hand point of the 
intersection with R. By the results of the preceeding paragraph it suffices to compute 
the probability that a Brownian motion started at f{xj + if e {xj)) will hit (p(L r ) before 
returning to T. By conformal invariance, this probability is comparable to f e {xf)/ f r (xj) 
a quantity which is equivalent by Koebe to 

r {(p'jxj + if r (x j )\ 
e \(p'(xj +if £ (xj)\' 

On the other hand, by quasiconformality, 

(x j+1 -x j )\tff{x j +if £ (x j )\ 

Combining all the estimates we see that the probability we look for is comparable to 



e d 



r\ x - lip'jxj + ifrjxj))] d (r^ ^ (x j+1 - Xj ) 



and we see a Riemann sum appearing: we finally get as an estimate for the probability 
we seek 

— length(<p(L r )) ~ - 

r r r a \Es 

which is precisely what we wanted. 



The preceeding proof has been presented because it is particurlarly simple, but of course 
the result is not optimal. The results of paragraph 2 remain true in dimension 2 and the 
proof requires only minor changes. In particular the result is true for all quasicircles; the 
difference with the case presented is that the topology of the level sets of the function 
distance to the boundary is more complicated in general. 




Figure 3: Two examples of self avoiding walk in 2D generated by the pivot algorithm 
fKenOS^ . The number of steps in each SAW is fixed at 10 5 . 

5 The self-avoiding walks case 

Self-avoiding walks (S.A.W.) (see |MS93] and [dG79] for a definition) serve as a good 
model for polymers in physics. On the other hand it is strongly believed that in 2D 
self-avoiding walks is the same as SLE 8 / 3 (see |RS05] for a definition of SLE K ). This 
conjecture, highly plausible, is comforted by the adequation between the computed 

4 

dimension, which is — and the proved dimension for SLE K curves , 1 + | (Beffara, 

|Bef06j ). The following simulations can be seen as a new way of probing the adequation 
between SAW's and SLE 8 / 3 : 

In order to check that the statistics of flights over a self avoiding walk follow the expected 
law with d — |, we have performed extended computer simulations. We have generated 
a set of self-avoiding walks on a square lattice using an implementation of the pivot 
algorithm described by Kennedy [Ken02j (see also [Ken05j ). The number of steps of 
the self avoiding walk is fixed at 10 5 . Two S.A.W. are shown in Fig. [3j 

We have checked by a box counting method that the mass fractal dimension of these 
curves are numerically founded around 1.33 ± 0.005. These values are very close to 
the expected value 4/3. We have performed an on-lattice simulation analyzing the first 
passage statistics of a random walk starting in the close vicinity of the SAW and going 
back for the first time nearby the SAW. The numerical computations were performed on 
several configurations of SAW, using a statistical analysis over more than 210 9 flights. 
Two probability density functions were computed. First, the probability density ip{n) 



that a flight has a total length equals to n. Second, the probability distribution of 
displacements 6(r). In order to limit edge effects, we have selected flights starting and 
ending on the same side of the S.A.W.. As shown Fig. [4] and Fig. [5l we found that 
ijj(n) oc n~ a and 8(r) oc r~@. 
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Figure 4: Evolution of the probability density 6(r) that a particle starting from close 
vicinity of SAW, returns to the SAW, for the first time, after an end to end displacement 
found between r and r + dr . The numerical estimation of the exponent p is very close to 
7/3. 



It was shown in |GKL + 06j that for a boundary of fractal dimension d embedded in 
an Euclidian space of Euclidean dimension d e , we should have 

d - d e + 4 

a = . (6) 

We are now in position to prove this estimate rigorously; this result is contained a paper 
to come. 

Moreover as 9(r) = —dP(r)/dr, we get 

p = d-d e + 3. (7) 
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Figure 5: Evolution of the probability density ip(n) that the particle starting from a 
close vicinity of SAW and returning for the first time to this SAW, has a total length 
displacement between n and n + 1. The numerical estimation of the exponant a is very 
close to 10/6. 

For d = 4/3, we expect to have a = 10/6 and (3 = 7/3. As shown in Figs. H^nd El 
numerical results provide a very good approximation of these above predictions. 

After these convincing simulations, let us prove rigourously these asymptotics for 
SLE K curves. First, combining results of Rohde-Schramm [RS05j and Beffara [Bcf06j, 
we see that the limsup in the definition of Minkowski upper- dimension for SLE K is actu- 
ally a limit, allowing asymptotic values for all values of r. Secondly, we observed in the 
last section that the result follows from the understanding of the number of Whitney 
cubes of given order. By a nice result of Bishop |Bis96] , it follows immediately that we 
get the right estimate if we allow the starting point to be chosen on both sides of the 
curve, which is actually the case in the above simulations. 

If we always start from the same side, then not only Bishop's result does not apply but 
neither does it follow from the previous discussion because the corresponding domains 
do not have the corkscrew property. But, as Rohde and Schramm have proved, these do- 
mains are Holder, meanning that the Riemman mapping from the upper half-plane onto 
them is Holder continuous. This condition implies a weaker form of the the corkscrew 



condition which is sufficient to ensure the possibility to compute the Minkowski dimen- 
sion of the SLE K curve via the counting of Whitney cubes. As we have seen this is 
enough to prove the main result about statistics of flights. 

The case of self affine curves is a little more delicate and will be treated in a forth- 
coming paper. 

It is to be noticed that, using quantum gravity arguments, Duplantier also obtained 
the right exponents for all SLE K (cf. |Dup04| ) 
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